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Abstract 

We use a simple system, the electron configuration in a Hydrogen-like atom, to demonstrate 
the importance of using a complete basis set to provide a proper quantum mechanical description. 
We first start with what might be considered a successful strategy — to diagonalize a truncated 
Hamiltonian matrix, written in a basis consisting of Hydrogen (Z = 1) basis states. This fails to 
provide the correct answer, and we then demonstrate that the continuum basis states provided the 
rest of the true wave function, for the bound ground states. This work then shows, in a relatively 
simple system, the need to utilize a complete basis set, consisting of both bound and continuum 
states. 
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I. INTRODUCTION 


Undergraduate texts on Quantum Mechanics generally emphasize analytical solutions 
to the Schrodinger Equation. But they almost always include a section on Hilbert Space, 
and basis states, and the formal solution of a problem through matrix mechanics. In a 
number of recent papers, ~ 1 we have illustrated, with familiar examples (Harmonic oscillator, 
Coulomb, etc.), the determination of bound state energies and eigenfunctions through a 
matrix formulation of quantum mechanics. All of these examples have used the familiar basis 
states which are the eigenstates for the infinite square well. In all the cases considered, we 
have truncated the Hilbert Space to a manageable size, so that the low-lying eigenvalues and 
their accompanying eigenfunctions are determined to high accuracy. 

At the same time, while not really emphasized in textbooks, students are generally made 
aware that a proper Hilbert Space needs to complete, to be useful as a basis set. This message 
was driven home recently in a study of the Helium electronic ground state, utilizing a basis set 
consisting of two-electron product states of Hydrogenic bound state eigenfunctions. While 
this is the natural basis set to use to understand Helium (and, indeed, any multi-electron atom 
of the periodic table), Hutchinson et al. emphasized that these product states constitute 
an incomplete basis set, and an accurate description requires the continuum states as well. 
Quantum chemists learned this lesson long ago, 1 and therefore never use such a (technically) 
poor basis set. As argued in Ref. [1], however, for state-of-the-art research problems on 
correlated electron systems, sometimes one does not have a choice. 

The case of the Helium electronic ground state is sufficiently complicated that the lesson 
in Ref. [ ] may be beyond the reach of undergraduates. Our purpose is to present a much 
simpler case, that of the one electron Hydrogen-like atom with central charge Zq, which we 
call ‘Z 0 Hydrogen’. Here of course the energies and states are a trivial extension to Hydrogen 
itself, where the Bohr radius a 0 is simply changed to a 0 /Z 0 . 

We start first by supposing that we do not know this answer, but instead decided to 
formulate the problem as a matrix problem, using as a basis the familiar bound states of the 
Hydrogen atom. There are an infinite number of these bound states, so we must necessarily 
truncate. As shown in the next section the result converges very rapidly, and no more than 
ten or so Hydrogen bound states are required to attain a converged result, but it is the wrong 
result ! The reason is explained, also in physical terms, and in the ensuing section we make 
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use of the known exact ground state for this problem to project out the contributions from 
each Hydrogenic state, including the continuum ones. These coefficients now sum to unity, 
as should be the case, for a proper description of Z 0 Hydrogen. We also illustrate how the 
continuum wave functions with increasing momentum contribute to the final (correct) wave 
function. 


II. ATTEMPT AT A MATRIX FORMULATION FOR Z 0 Hydrogen 


We are interested in the ground state solution for the problem of a single electron 
interacting with a fixed nucleus of charge Z 0 . This has a Hamiltonian given by 

e 2 


H = h 


V 2 -Z n 


.2 1 


( 1 ) 


2m 47reo r 

where U is Planck’s constant, m is the mass of the electron, —e is the charge of the electron, 
eZ 0 is the charge of the nucleus, eo is the permittivity of free space, and r = |rj is the 
radial coordinate. The exact eigensolutions to this problem acquire the usual three quantum 
numbers, (n, £, m), and as for the case of Hydrogen, the ground state has n — 1 and t = m = 0, 
and has energy = —ZqE[°\ where E^ = fr 2 /(2mal) = 13.6eV, and a 0 = ine 0 h 2 /(me 2 ) 
is the Bohr radius. We have assumed the nucleus to be infinitely heavy so the same mass 
appears in Eq. (1) as appears in the definition of the Bohr radius. The exact ground state is 
given by' 1 ’ 


= -i= V e~ Zor/ao 
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where the is the contribution from the angular part of the wave function for all 

s-states. 

Now if we proceed as described in the Introduction, pretending not to have knowledge of 
this result, we would first expand the (unknown) ground state wave function in terms of a 
‘handy’ set of basis states, which we denote as </h, with the index in principle representing 
multiple quantum numbers. Since we are using a central potential, and we anticipate the 
ground state to have s-wave symmetry, then all the basis states have this symmetry as well. 
Thus, if we use the Hydrogen bound states as a basis set, then the label l n , will denote the 
principal quantum number n, and £ — m — 0. Writing 


'y ^ i0 n ), 
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and taking inner products with each (orthonormal) basis state and the Schrodinger Equation 
written in this basis results in the matrix equation 

OO 

^ ) HnmCm Ec m . (4) 

771=1 

This represents an infinite dimensional matrix equation, so to make progress we truncate at 
n m axi vary this maximum number, and monitor the convergence of the ground state energy, 
for example. The required matrix elements are 


Hnm {(pn\H\(j) m 'j. 


( 5 ) 


The simplest way to proceed is to rewrite the Hamiltonian, Eq. (1), as 
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where H 0 is the actual Hamiltonian for Hydrogen (first two terms in second line) and H’ 
is the remaining term. We begin with the diagonal terms, (0 n |f/|0 n ). Because 0 n is an 
eigenstate of the Hydrogen hamiltonian, these terms simplify drastically: 
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Here we have used the results for the energy levels of hydrogen (E n = —and the well 
known result (0„|7|0 n ) = 1 /(n 2 a 0 ). 

Because the individual 0 n are orthonormal eigenstates of H 0 , the off-diagonal terms reduce 
to a simple inner product, 


(0n|7fo + 7^1 0m) — E n {(j) n \ 0 m ) + (0n| 7T|0m) 

= -{Zo - dj— (0n|-|0m) (8) 
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= —2(Z 0 — l)E’| O) (0 n | —|0 m ). 
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Knowing the eigenstates for Hydrogen, these integrals are individually straightforward. A 
general formulation requires the wavefunctions of hydrogen: 
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where in the second line l = m = 0. These are then substituted into the inner product in 
Eq. (8). The angular integral can be done immediately, as there is no angular dependence, 
and this eliminates the -^=; furthermore, note that hereafter the letter ‘m’ denotes a principal 
quantum number (not the azimuthal quantum number which is now always zero): 
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where I nm is simply a number. This integral can be done numerically. Alternatively, an ana¬ 
lytic solution to the integral is achieved by writing out the associated Laguerre polynomials 
as a (finite) power series. Using this series makes the required integral elementary, so we 
end up with 
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which is simply a number. Thus, the off-diagonal results for the Hamiltonian are: 


(<f>m\H\<f> n ) = -2 (Z - l)E} O) (0 m | —|0 n ). (12) 

r 

Following the philosophy of Ref. [i], we can simply determine the Hamiltonian matrix 
up to some maximum cutoff, n max , and diagonalize it to determine the ground state energy. 
The results are shown in Fig. 1. In Fig. (la) we see that as we increase the size of the 
matrix for different values of Zq, the ground state energy does converge (almost immediately). 
However, the energies converge to the wrong value (except the case of Z 0 — 1, which is 
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simply hydrogen, and obviously needs only the one basis state for the correct answer). The 
expected value is also indicated - it is simply —ZqE[°\ as noted previously. In Fig. (lb) we 
show the actual energy achieved vs. Zq, along with the exact result, and their difference. 
These clearly diverge, especially as Zq increases. 
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FIG. 1: (color online) (a) Ground state energy (in units of Rydbergs) for various values of Zq 
as a function of matrix size. Note that in all cases convergence is achieved with n max ~ 5. (b) 
Ground state energy (in units of Rydbergs) vs. Zq (solid blue curve) using n max = 18. Also shown 
is the exact result, Eq/e[°^ = — Zq (green dashed curve) and the difference (dotted red curve). 
The deviation grows as Zq increases, but notice some small discrepancies for Zq < 1 as well, (c) 
An expanded view of the Zq = 0.5 result from part (a), to better illustrate the discrepancy with 
the exact result, (d) An expanded view of the dependence on Zq, for Zq < 1 , showing that the 
difference peaks near Zq ~ 0.6. Note that 10 x the difference is plotted to see it better. 


What went wrong? As already mentioned in the Introduction, the bound state eigenstates 
for Hydrogen, while infinite in number, do not actually form a complete basis set. We should 
not have expected to get the correct result. On the other hand, truncating the Hilbert space 
has worked in previous studies of one and three dimensional problems. 1,3 The problem here 
is that as we increase Zq, we are trying to describe an atom whose electron is more tightly 
bound than in Hydrogen (note that (r) oc Oq/Zq). But utilizing the bound excited states 
attempts to make use of states that are more extended, not less extended. Note that for 
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Zo < 1 the calculation is more accurate. However, in Fig. (lc) we show an expanded view of 
the result for Z 0 = 0.5, and even here the matrix result disagrees with the exact result by 
a small but definite amount. A closer examination of the entire low Z 0 region is provided 
in Fig. (Id) and we see that the error has a maximum about halfway between Zq = 1 and 
Zq = 0. The relative success in this region is because physically we are trying to construct 
a state that is more extended compared to Hydrogen, and the excited states are helpful in 
producing this. In contrast, for Z 0 > 1 we require the continuum states as well, as they are 
perfectly capable of describing a more closely bound state (recall that an infinite set of plane 
waves can describe a h-function). We demonstrate this in the next section. 

It should be noted that for the sake of pedagogy a different tact could have been followed. 
One could imagine wanting to solve for the ground state of ZoHydrogen, not in terms of 
the bound eigenstates of Hydrogen, but in terms of Hydrogenic states with nuclear charge 
Z re f, where the subscript ‘ref’ is short for ‘reference’. This is more artificial then the problem 
already considered, but it is instructive to consider anyways, and we do so in the Appendix. 

III. THE CONTRIBUTION FROM THE CONTINUUM AND A FULL SPEC¬ 
TRAL DECOMPOSITION 

Once we recognize the need to include the continuum states in the correct description, it 
immediately becomes difficult to formulate this problem as a finite-sized matrix diagonaliza- 
tion problem. Here, however, we know the exact solution, so it is still possible to determine 
the degree to which each basis state contributes to the overall solution; this will clearly vary 
with Z 0 . With the radial part of the wave function denoted by (r) = V^V ; K) 0 o( r ) [ see Eq. 
(2)], the correct spectral decomposition is given by 

OO 

r) = ^a n R n o(r ) + 

n= 1 

where the a n and a v coefficients can be determined by overlap integrals (assuming the 
left-hand-side is known), and these refer to the discrete (n) and continuum (p) components, 
respectively. Note that in both cases symmetry considerations require only l — 0. 

The procedure will be most familiar for the discrete coefficients, so we begin with these. 
We take overlap integrals of Hydrogenic wave functions (with Z — 1) with the exact wave 
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function, i.e. 
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The radial wave function R n o(r) can be written in terms of an associated Laguerre polynomial, 
which has a known finite polynomial expansion; then the integral in Eq. (14) can be readily 
done, and what remains are finite sums which one can recognize as binomial expansions. The 
final result is 
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Note the explicit factor of (Z o — 1) in front of the second term; for Zq = 1, the only non-zero 
coefficient is the n — 1 term, with coefficient unity, as must be the case. For Z 0 ^ 1, all other 
s-states contribute with diminishing amplitude as n increases. 

A numerical summation of the probabilities \a n \ 2 over all values of n reveals that these do 
not sum to unity (when Zq ^ 1), i.e. a finite contribution must come from the continuum 
states. The continuum states for l = 0 are denoted R p o(r), where p is a continuum momentum. 
These eigenstates for the Hydrogen atom are less familiar to students, but they are given in 
standard undergraduate texts: 1 ’ ‘ 
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where we have written this for general Z but require only Z — 1, and 


M(a,b,z ) = 
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is the so-called Rummer function, ■ and (a) m = a(a + l)(a + 2)...(a + m — 1) is the so-called 
Poc hh a mm er symbol. Note that (a)o = 1, and (2) m = (m + 1)!, for example. Equation (16) 
is a general solution, and is connected to the Laguerre polynomials that describe the radial 
bound state wave functions. 

The standard normalization condition for the continuum states, 


dr r 2 Rpi(r)R p (r) = 5(p' — p), (18) 

determines the coefficient in Eq. (16). By utilizing the expansion in Eq. (17) we can evaluate 
the overlap integral required to obtain the coefficient a p , for Z — 1, 



dp 


4 (Z 0 a 0 ) 3/2 (Z 0 - 1) / 2-iTpao 

[(pa„r + ZSP Y ! _ e -* exp 



(19) 














0.6 



FIG. 2: (color online) The momentum distribution a p /ao given in Eq. (19) as a function of pa® for 
various values of Zq. For Zq < 1 there is a negative correction, i.e. the continuum states try to 
make the approximation based on the bound states more extended, while for Zq > 1 the corrections 
are all positive, so try to increase the value of the wave function near the origin. The scale on which 
the continuum momentum eigenstates contribute most increases with increasing Zq. 

These coefficients, when squared, show a distribution peaked around p ~ Z 0 /a 0 , as expected. 
Since these enter linearly in the wave function expansion given in Eq. (13), we show in Fig. (2) 
the amplitudes, a p /ao vs. pao for a variety of values of Zq, showing how contributions from 
the continuum peak near pao ~ Zq. The larger Zq, the larger is the contribution from the 
continuum states. 

To sum the continuum contributions one requires a factor of 2/tt to account for the 
enumeration of continuum states; then the contribution from the all the continuum states is 

p “"* = ^ 327 , ~zi~ l y) ' ( 0) 

The contributions from the various states to the ground state are displayed in Fig. (3). 
We have checked that for all Zq the contributions sum to unity. Note that for Zq > 1 the 
continuum states play an increasingly important role. Let us repeat the reason, now that we 
Fig. (3). In this regime the actual wave function, given by Eq. (2), varies on a scale of ao/Zo. 
For Z 0 > 1 none of the bound states can provide structure on such a scale; as n increases the 
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Z 0 

FIG. 3: (color online) Various contributions from the Hydrogenic bound and continuum states to 
the ground state wave function for a central charge with magnitude eZo, as a function of Zq. At 
Zq = 1 only the Is state contributes, as expected. For Zq ^ 0 all other states contribute as well. 
In particular, for Zq > 1 the continuum states contribute with increasing amplitude for reasons 
explained in the text. The sum of the three curves is unity for all Zq, as indicated by the thick 
horizontal (purple) line across the top. 


Hydrogenic basis states become more extended, not less. The only source of variation on 
this scale is the spectrum of continuum states, and these are readily utilized. 

The degree to which different continuum momentum eigenstates contribute is illustrated 
in Fig. (4), where we show the wave function given by Eq. (13) for a particular example, 
Zq = 1.5, but with the momentum integration cut off at increasing values of momentum. As 
compared with the exact wave function it is clear that initial components (near pciQ ~ 0) first 
fix the large r behavior (the bound states produced a value that was too high for r/ao > 1) 
and components with larger value of momentum then begin to produce larger amplitude 
near the origin. This figure demonstrates explicitly what we were just explaining, that finer 
scale contributions were necessarily produced by the continuum states. 

This example serves to show how a ‘poor’ basis choice (i.e. using the Hydrogenic set for 
Z = 1 rather than for Z = Zq) can require the use of the continuum states. For the case of a 
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r/a 0 

FIG. 4: (color online) The evolution of the ground state wave function from just considering the 
bound states, i.e. the lowest (blue) curve to the exact result (blue squares). As the accuracy 
increases we cut off the momentum integration in Eq. (13) at p c ao = 1 (next highest (red) curve), 
PcCio = 3 (next highest (green) curve), p c ao = 10 (next highest (black) curve), and finally, p c ao = 
50 (next highest (red) curve). The higher momentum continuum components contribute to the 
wave function near the origin. 


single electron bound to a positive charge, a ‘poor’ choice can obviously be avoided. However, 
in the case of many electron problems where locally one can have either a single electron or 
two electrons, it is very difficult to devise a basis set that diagonalizes both scenarios. 1 


IV. SUMMARY 

We have used one of the simplest yet realistic quantum systems, the electron configuration 
in a Hydrogen-like atom, to demonstrate a principle with which most students are familiar, 
but likely few have encountered in practice: the need for a complete basis set to describe 
properly a quantum mechanical system. The Hydrogen-like atom serves as a good system to 
show this, since we know the exact solution, and hence know in advance the correct answers. 
We can also do most of the required integrals, both for the matrix formulation, and for the 
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spectral decomposition. The matrix formulation did not work, perhaps counter to what 
some might expect. The implication of this failure, that continuum states are necessary, 
was reinforced by carrying out a spectal decomposition. This result, portrayed in Fig. (3), 
perhaps runs counter to our intuition, that bound states require a partial amplitude (that 
can be substantial!) corresponding to continuum states, i.e. states in which the particle is 
free to roam through all space. We also provided a physical understanding of why these 
states were especially required when we attempt to construct states that are more tightly 
bound than provided by the bound basis functions, as demonstrated explicitly in Fig. (4). 
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Appendix A: A Matrix formulation using Z ve f 7 ^ 1 


As in Section [II] we wish to solve for the Hamiltonian 

— fi 2 2 e 2 1 
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but we now rewrite this using a slightly different decomposition: 
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The motivation is that we need more tightly bound wave functions to describe ultimately 
the ground state for Z 0 Hydrogen. But if Z ref » Z 0 then many of the bound states will 
actually be more ‘compact’ compared to the ground state we are seeking to find. Therefore, 
proceeding as before, we begin with the diagonal terms, {(p n \H\(j) n ), to obtain 
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Here we have used the results for the energy levels of a hydrogenic state with nuclear charge 

E ( 0 ) 

Z iei : E n = — Z 2 e{ -^-, and the generalization of the well known result 1, {<f> n \b\<f>n) — Z ie {/(n 2 a 0 ). 



FIG. 5: (color online) The ground state energy for a ZoHydrogen system with Z$ = 3, vs. Z re f. 
This calculation was performed by diagonalizing a 18 x 18 matrix, and shows that using larger 
values of Z re f is effective to a degree. Results are definitely more accurate than using Z re f = 1 for 
example, but the convergence to very accurate results is very slow. 


Evaluation of the off-diagonal matrix elements proceeds as before; the Enal result is 
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is the same number given by Eq. (10) in Section II. 

Now we can simply diagonalize an N x N matrix for some Z ve f. In Fig. (5) we show the 
ground state energy obtained in this way for a targeted system with Zq = 3. Use of the usual 
Hydrogen states with Z ref = 1 produces a very poor result, E 0 ~ —5 E^\ compared to the 
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exact result, E$ x = —9 E[ [) \ As we increase Z Te f the result becomes more accurate, and of 
course exact as Z ve f —> Z 0 . For larger values of Z ref the result hrst deteriorates before steadily 
improving as Z re f continues to increase towards 100, as shown. While the convergence is 
quite slow, nonetheless this exercise demonstrates that using large values of Z re f allows one 
to get reasonably accurate results for the reason that the bound state basis set now contains 
quite a number of wave functions that describe a particle confined to the origin. Nonetheless, 
the lack of better accuracy is an indication that a continuous set of continuum states in the 
end offers more flexibility than a discrete set of bound states. Note that we did not require 
values of N in excess of about 20 to achieve converged results for all values of Z re f shown. 
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